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1.0  ABSTRACT 


A  need  exists  for  optical  target  recognition  in  the  presence  of  clutter.  Our  innovation  has 
been  demonstrated  to  improve  the  signal-to-background  ratio  of  the  target  correlation  spike  by 
masking  off  most  of  the  cluttered  input  scene  to  reduce  the  background  noise.  The  correlation 
process  is  repeated  for  each  small  section  of  the  input  scene  to  find  the  visible  pieces  of  the  target, 
and  is  then  repeated  one  final  time  with  the  mask  passing  all  the  sections  of  the  input  scene  which 
were  found  to  have  pieces  of  the  target  but  without  the  surrounding  clutter. 

The  significance  here  is  found  in  the  fact  that  for  heavily  cluttered  targets  a  higher 
correlation  spike-to-background  ratio  is  obtained  by  reducing  the  size  of  the  input  scene,  even  if 
visible  parts  of  the  target  are  left  out.  This  occurs  because  the  input  scene  is  mostly  background 
clutter  with  a  few  separated  pieces  of  visible  target.  This  Consensus  Filter  mimics  the  human 
capacity  to  detect  targets  in  heavily  cluttered  scenes  by  concentrating  on  small  sections  and  ignoring 
the  surrounding  areas.  It  also  overcomes  the  poor  repeatability  of  a  human  who  is  working  on 
large  amounts  of  input  data.  The  new  requirement  of  having  to  repeat  the  correlation  for  each  piece 
of  the  input  scene  simply  better  utilizes  the  extreme  speed  of  the  optical  correlator.  The  tradeoff 
that  this  new  technique  entails  is  the  need  to  evaluate  the  large  number  of  spikes,  each  of  which 
indicates  a  possible  piece  of  the  target. 

The  demonstration  involves  a  military  scout  car  as  a  target  shape.  The  target  is  mostly 
obscured,  leaving  only  small  portions  visible.  The  obscured  areas  are  replaced  by  noise  and  the 
Consensus  Filter  is  applied  to  the  whole  input  scene.  The  computer  model  of  the  optical  system 
breaks  the  input  scene  into  blocks  of  4  x  4  pixels,  8x8,  16  x  16,  and  32  x  32  pixels.  We 
observed  that  the  Consensus  Filter  generally  has  a  higher  probability  of  detection  than  a  standard 
correlator  for  highly  obscured  targets.  In  addition  the  ratio  of  the  correlation  spike  to  the  next 
largest  spike  is  always  larger  using  the  consensus  filter  because  of  the  reduction  in  noise 
background  and  in  some  cases  by  identifying  a  more  optimum  fourier  plane  filter. 


Further  results  for  phase-only  filters  includes  the  effects  of  continuously  changing  the 
amplitude  information  in  the  fourier  plane  filter  from  the  matched  filter  condition  to  the  phase-only 
condition,  and  beyond  to  a  “hyper”  phase-only  filter.  This  sequence  continuously  increases  the 
performance  in  terms  of  the  ratio  of  the  main  correlation  spike  to  the  secondary  spike. 

In  addition  a  study  was  undertaken  to  find  an  acceptable  quality  measure  of  a  correlation 
plane  to  be  used  internally  by  the  Consensus  Filter.  One  measure  is  (max-mean)/mean.  The  mean 
and  standard  deviation  of  this  measure  were  studied  as  a  function  of  the  target  size  that  was  used  to 
create  the  phase-only  filters,  and  this  measure  was  successfully  utilized  in  the  Consensus  Filter. 
Further  work  is  needed  on  a  figure  of  merit  to  make  the  Consensus  Filter  a  fully  self  contained 
alternative  to  the  standard  correlator. 


Other  results  of  this  research  include  comparisons  of  matched  filters  and  phase-only  filters _ 

for  partially  obscured  targets,  and  for  random  target  shap«s.  .,..-,uiaauy  codaa 
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3.0  INTRODUCTION 


There  is  a  need  for  target  recognition  in  the  presence  of  clutter.  Optical  correlation  has  the 
advantage  of  achieving  great  speed  through  the  parallel  nature  of  the  optical  fourier  transform,  and 
is  well  adapted  to  the  Consensus  Filter  which  requires  a  large  number  of  optical  correlations. 

Our  version  of  the  clutter  problem  is  attempting  to  recognize  a  target  which  is  partially 
obscured,  say  by  a  stand  of  trees,  but  which  is  otherwise  undistorted. 

3.1  Description  of  Results 

The  demonstration  of  the  Consensus  Filter  was  accomplished  by  computer  simulation  of 
the  optical  correlator,  including  the  obscured  input  target,  the  fourier  plane  filter,  and  the  ancillary 
processing  based  on  data  taken  from  the  correlation  plane.  The  software  programs  were  written  by 
the  principal  investigator  using  Matlab  which  is  a  matrix  manipulation  language,  and  were  run  on 
personal  computers  using  8086,  80386,  and  80486  microprocessors. 

The  results  of  this  work  are  that  the  Consensus  Filter  concept  does  better  succeed  at  finding 
pieces  of  the  target  than  an  ordinary  correlator  when  most  of  the  target  is  missing.  An  ordinary 
correlator  is  superior  when  most  or  all  of  the  target  image  is  available.  In  addition  the  Consensus 
Filter  is  significantly  more  complex  and  will  require  additional  development  to  optimize  its  internal 
decision  making  for  a  wider  class  of  target,  clutter,  and  noise  scenarios. 

An  essential  component  to  the  operation  of  the  Consensus  Filter  which  involves  comparing 
target  correlation  spikes  from  different  output  correlation  planes  is  a  universal  figure  of  merit  This 
issue  was  investigated  and  partially  solved,  permitting  a  proper  demonstration  of  the  Consensus 
Filter.  This  work  is  discussed  in  Section  5. 

The  general  issues  of  phase-only  or  partially  phase-only  filters  were  partially  investigated 
with  a  focus  on  the  Consensus  Filter  application.  An  interesting  viewpoint  emerged  from  a  series 
of  tests  reported  in  Section  6  which  shows  that  phase-only  filter  acts  more  strongly  than  a  matched 
filter  to  spectrally  “whiten”  the  input  target  and  therefore  produce  a  correlation  spike  that  better 
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resembles  an  impulse.  This  target  whitening  is  extended  using  a  “hyper”  phase-only  filter  which 
does  have  amplitude  components  and  produces  a  correlation  spike  which  is  the  closest  match 
possible  to  a  perfect  impulse. 

3.2  Review  of  the  Consensus  Correlation  Concept 

The  Consensus  Correlator  operates  by  chopping  the  input  scene  into  pieces  which  are 
individually  checked  by  a  standard  optical  correlator  which  is  looking  for  the  whole  target.  The 
blocks  of  the  input  scene  which  are  found  to  contain  pieces  of  the  target  are  checked  to  be  sure  that 
their  correlation  spikes  are  close  together,  which  ensures  that  the  various  target  pieces  that  have 
been  found  would  fit  together  properly. 

The  essential  element  here  is  that  the  size  of  the  chopped  up  blocks  from  the  input  scene 
which  are  individually  handled  by  the  optical  correlator  must  be  smaller  than  the  perfect  target 
shape  which  was  used  to  create  the  fourier  plane  filter.  This  condition  results  in  the  background 
level  of  the  correlation  plane  being  reduced  in  the  vicinity  of  the  correlation  spike.  The  Consensus 
Correlator  is  based  on  the  idea  that  only  a  small  portion  of  the  target  is  present  and  that  cutting  back 
on  the  noise  or  clutter  near  this  small  target  piece  will  give  a  cleaner  correlation  spike.  The 
relevance  of  the  target  size  can  be  understood  from  the  convolution  viewpoint  where  noise  and 
clutter  are  spread  a  distance  in  the  correlation  plane  equal  to  the  width  and  height  of  the  target 
shape.  (This  is  approximately  true  also  for  typical  phase-only  filters.)  For  example  consider  the 
input  scene  cut  into  blocks  which  are  10()  times  smaller  in  area  than  that  of  the  target.  The 
correlation  plane  receives  the  uncorrelated  (noise)  optical  power  spread  out  by  a  factor  of  100.  If 
the  correlation  spike  is  unchanged,  the  signal  to  background  ratio  in  the  vicinity  of  the  correlation 
spike  will  be  100  times  greater. 

The  improvement  in  the  signal  to  background  ratio  occurs  because  the  fourier  plane  filter 
spreads  out  the  noise  or  clutter  from  the  small  piece  of  the  input  scene.  If  the  masked  off  block  of 
the  input  scene  is  too  small  then  the  correlation  spike  will  also  be  reduced,  because  the  small  piece 
of  the  target  which  is  visible  in  the  input  will  also  be  cut  up. 


4 


The  preference  for  cutting  the  input  scene  into  blocks  roughly  equal  to  the  size  of  the 
expected  target  pieces  is  accommodated  by  running  the  system  with  several  different  block  sizes 
and  choosing  the  strongest  performer.  Our  demonstration  used  4x4  pixels,  8x8,  16  x  16,  and 
32  x  32  pixels.  We  observe,  as  expected,  better  performance  when  the  input  is  chopped  into 
blocks  closest  to  the  size  of  the  target  pieces  used  in  the  input  scene  to  the  correlator. 
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4.0  CONSENSUS  CORRELATOR  DEMONSTRATION 


We  have  successfully  demonstrated  the  Optical  Consensus  Correlator  by  computer 
simulation.  This  demonstration  consists  of  properly  locating  a  highly  obscured  target  with  a  higher 
probability  than  is  obtained  by  using  a  traditional  optical  correlator.  The  origin  of  the  target  shape 
is  a  military  scout  car  shown  in  Figure  1 .  The  target  is  represented  as  zeros  and  ones  on  a  32  x  32 
matrix,  as  shown  in  Figure  2.  The  actual  input  image  used  is  produced  by  masking  off  all  but  a 
few  pieces  of  the  target  matrix,  and  replacing  the  masked  areas  by  noise,  which  represents  clutter. 

The  noisy  pixels  are  statistically  independent  and  are  generated  as  the  absolute  value  of  a 
Guassian  distributed  zero  mean  process  with  standard  deviation  equal  to  2  for  the  maximum  clutter 
level  considered.  Note  that  input  to  the  correlator  is  in  terms  of  electric  field  and  the  correlation 
plane  is  dealt  with  in  optical  power.  Finally,  the  graphs  showing  the  probability  of  detection 
versus  noise  level  range  in  noise  level  from  0  to  10;  10  representing  a  noise  standard  deviation  of  2 
for  the  electric  field  input,  and  a  1  representing  a  noise  standard  deviation  of  0.2  for  the  electric 
field  input.  The  optical  noise  or  clutter  power  increases  as  the  square  of  the  number  given  on  the 
x-axes. 


The  Consensus  Correlator  uses  a  scanning  mask  to  restrict  the  area  being  evaluated  in  the 
input  scene.  Four  different  mask  aperture  sizes  are  used,  and  the  target  detection  probability  is 
given  for  each  one  in  the  graphical  data:  4x4  pixels  blocks,  8x8,  16  x  16,  and  32  x  32  which 
gives  the  response  of  the  standard  correlator  for  comparison. 

4.1  The  Consensus  Correlator  Algorithm 

This  section  is  a  block  description  of  the  algorithm  used  to  demonstrate  the  Consensus 
Filter,  which  is  reproduced  in  full  in  the  Appendix. 

4.1.1  Scanning  the  Input  Scene:  scan2.m 

The  input  scene  is  a  matrix  which  is  first  masked  off  except  for  a  small  block  as  shown  in 
Figure  3,  step  1.  The  masked  off  input  matrix  is  fourier  transformed  and  then  (array)  multiplied  by 
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CONSENSUS  FILTER 


the  phase-only  target  filter,  inverse  transformed,  absolute  valued  and  squared  to  produce  the 
correlation  plane.  The  largest  spike  is  noted  and  other  spikes  are  noted  if  they  are  above  a  certain 
threshold  fraction  of  the  largest  spike  (software  variable  :  ff).  The  maximum  number  of  spikes 
accepted  is  limited  (variable  :  spkbx).  The  location  of  these  spikes  is  recorded.  This  analysis  is 
repeated  for  the  entire  input  scene  using  one  piece  of  the  input  scene  at  a  time. 

4.1.2  Grouping  the  Correlation  Spikes:  test3.m 

The  purpose  here  is  to  determine  if  there  is  any  overlap  of  the  correlation  spikes  from  the 
various  sections  of  the  input  scene.  An  overlap  suggests  that  the  correlation  spikes  were  produced 
by  various  parts  of  a  target  that  would  fit  together  properly.  The  algorithm  hunts  for  closely 
located  correlation  spikes,  with  an  adjustable  maximum  separation  (variable  :  space). 

4.1.3  Final  Mask  Comparison:  finmask2.m 

This  step  creates  a  new  mask  which  is  placed  in  front  of  the  input  scene  and  passes  only 
those  blocks  of  the  input  scene  which  had  individually  produced  a  correlation  spike  at  a  common 
location.  See  Figure  3,  step  2.  This  masked  input  is  fed  to  the  correlator  and  the  resulting 
correlation  spike  is  located  and  measured.  Running  the  correlator  on  an  input  scene  where  only 
those  parts  of  the  scene  that  have  pieces  of  the  target  are  passed  will  result  in  a  much  larger 
correlation  spike  and  a  minimum  of  noise  or  clutter. 

A  further  improvement  occurs  at  this  point.  A  phase-only  filter  is  not  optimum  for  a 
partially  obscured  target;  and  this  is  discussed  at  greater  length  in  Section  6.  A  new  phase-only 
fourier  plane  filter  is  produced  based  on  the  apparent  location  of  the  target  (correlation  spike)  and 
the  implied  parts  of  the  target  found  in  the  input  scene.  (This  uses  optfilt2.m.)  This  new 
optimized  phase-only  filter  is  applied  to  the  masked  input  and  the  new  correlation  spike  observed. 
This  optimization  procedure  is  then  repeated  a  second  time.  This  phase-only  fourier  plane  filter 
which  is  optimized  to  match  the  visible  pieces  of  the  target  can  noticeably  sharpen  the  resulting 
correlation  spike. 

During  the  process  of  reconstructing  a  more  closely  matched  fourier  plane  filter,  the  mask 
is  also  checked  to  be  sure  it  doesn’t  pass  any  parts  of  the  input  scene  where,  based  on  the  new 
knowledge  of  the  target  location,  there  could  not  be  any  target  components. 
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This  whole  process  is  repeated  for  the  other  likely  mask  combinations  which  were  indicated 
by  groups  of  overlapping  correlation  spikes  from  the  first  scan  of  the  input  scene.  The  single 
optimized  correlation  spikes  resulting  from  the  resulting  mask  combinations  are  compared  and  one 
target  location  is  selected.  The  figure  of  merit  used  to  make  this  comparison  is  discussed  in 
Section  5. 

The  algorithm  is  designed  to  reduce  transmission  of  clutter  from  the  input  scene  to  the 
immediate  vicinity  of  the  detection  spike  in  the  correlation  plane.  This  concept  is  also  presented  by 
Figures  3  and  4. 

4.2  Four  Examples  of  Cluttered  Target  Detection 

4.2.1  Two  Small  Pieces  of  Target  Visible 

Figure  5  shows  a  contour  of  the  target  shape  with  two  4x4  pixel  boxes  drawn  in  to 
indicate  the  two  small  pieces  of  the  target  that  will  be  visible  at  the  input  of  the  correlator.  The 
remainder  of  the  matrix  has  noise  added.  Figure  6  compares  the  probability  of  detection  of  the 
standard  correlator  (solid  line)  with  the  Consensus  Correlator  (4x4  blocks  are  dash-dot,  8  x  8  are 
dot-dot.  16  x  16  are  dash-dash)  as  a  function  of  increasing  background  noise.  The  performance  of 
the  4  x  4  pixel  mask  used  to  scan  the  input  scene  gives  the  best  performance  as  is  expected  since  it 
is  a  good  match  with  the  size  of  the  visible  pieces  of  the  target.  The  probability  of  detection  is 
obtained  by  running  ten  different  versions  of  the  noise. 

Figure  7  shows  the  average  ratio  of  the  correlation  spike  to  the  secondary  peak  for  the  four 
cases,  with  zeros  averaged  in  if  the  target  was  not  correctly  located.  The  Consensus  Correlator 
both  reduces  the  noise  and  clutter  in  the  correlation  plane  and  optimizes  the  detection  filter,  which 
is  shown  to  improve  the  ratio  of  the  correlation  spike  to  the  secondary  peak:  ratio  = 
(max-median)/(2nd  max-median). 

4.2.2  One  Small  Piece  of  Target  Visible 

Figure  8  presents  the  target  with  an  emphasized  box  which  shows  the  only  visible  part  of 
the  target  with  noise  everywhere  else  in  the  input  matrix.  Figure  9  shows  the  probability  of 
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AREAS  OF  CLUTTER  IN  INPUT  SCENE 
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Figure  5.  A  contour  of  the  target  shape  with  two  4x4  pixel  boxes  to  emphasize  the  two  small 
pieces  of  the  target  that  are  presented  to  the  input  of  the  Consensus  Correlator. 


Figure  6.  The  probability  of  target  detection  shows  the  expected  decrease  with  increasing 

noise  level.  The  background  optical  noise  rises  as  the  square  of  the  x-axis  values  of 
0-10.  The  solid  line  is  the  standard  correlator.  The  scanning  16  x  16  pixel  box  is 
dash-dash,  the  8  x  8  is  dot-dot,  and  the  4  x  4  is  dash-dot.  The  4x4  mask  shows 
performance  superior  to  the  standard  correlator.  This  is  based  on  ten  runs  using 
different  noise  inputs. 
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Figure  7.  The  ratio  of  the  amplitude  of  the  primary  correlation  spike  to  the  amplitude  of  the 
secondary  peak  shows  an  advantage  of  the  Consensus  Correlator  in  reducing 
background  noise.  The  dash-dot  curve  is  the  4  x  4  pixel  mask  and  the  solid  curve 
is  the  standard  correlator.  The  horizontal  axis  is  the  background  noise  level. 


detection  as  a  function  of  noise  and  clutter  with  the  4  x  4  mask  (dash-dot)  giving  the  best  result. 
Again,  a  4  x  4  mask  is  a  good  match  with  the  size  of  the  visible  target  piece.  The  16x16  mask 
system  failed  even  at  low  noise  which,  due  to  a  few  similar  failures  elsewhere,  is  assumed  to  be  a 
software  problem. 

4.2.3  One  8x8  Pixel  Piece  of  the  Target  Visible 

Figure  10  shows  the  target  shape  with  an  8  x  8  pixel  box  emphasizing  the  small  fraction  of 
the  target  shape  which  is  visible.  The  remaining  area  of  the  input  scene  is  filled  with  noise. 

Figure  1 1  shows  the  probability  of  detection  as  a  function  of  noise  and  clutter  intensity.  The  8x8 
pixel  mask  (dot-dot)  gives  the  best  performance,  and  it  is  the  best  match  for  the  size  of  the  visible 
target  piece. 
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Figure  8.  The  emphasized  4x4  pixel  box  shows  the  smaller  portion  of  the  target  shape  used 
as  the  input  for  Section  4.2.2. 
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Figure  9.  Probability  of  target  detection  versus  noise  level  for  the  previous  figure.  The  4x4 
pixel  scanning  mask  (dash-dot)  is  superior  to  the  standard  correlator. 


Figure  10.  The  emphasized  8x8  pixel  box  shows  the  portion  of  the  target  shape  used  as  the 
input  for  Section  4.2.3. 
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Figure  1 1 .  Probability  of  target  detection  versus  noise  level  for  the  previous  figure.  The  8x8 
pixel  scanning  mask  (dot-dot)  is  superior  to  the  standard  correlator. 
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4.2.4  One  8x8  Pixel  Piece  of  the  Target  Visible  with  a  Shift  in  the  Target  Position 

Figure  12  shows  the  target  shape  with  an  8  x  8  pixel  box  denoting  the  fraction  of  the  target 
shape  that  is  visible.  In  addition  a  shift  in  the  position  of  the  target  was  arranged  to  highlight  the 
sensitivity  of  the  Consensus  Correlator  to  fortuitous  overlaps  of  the  scanning  mask  and  any  visible 
piece  of  the  target.  The  algorithm  has  not  yet  been  modified  to  desensitize  it  to  this  condition. 
Figure  13  shows  the  probability  of  detection  versus  noise  and  clutter.  In  this  case  the  4  x  4  pixel 
mask  (dot-dash)  is  superior  to  the  8  x  8  mask  (dot-dot)  despite  the  8  x  8’s  better  size  match.  The 
8x8  mask  apparently  accepts  excessive  noise  in  comparison  to  the  4  x  4  mask  due  to  a  weak 
overlap  in  the  mask  scanning  step. 

The  Consensus  Correlator  demonstrates  a  superiority  over  the  standard  correlator  when  a 
small  fraction  of  the  expected  target  is  visible. 
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Figure  1 2.  The  emphasized  8x8  pixel  box  shows  the  portion  of  the  spatially  shifted  target 
shape  used  as  the  input  for  Section  4.2.4. 


Background  Noise 

Figure  1 3.  Probability  of  target  detection  versus  noise.  The  shifted  input  target  from  the  Figure 
above  has  interfered  with  the  fortuitous  overlap  of  the  scanning  8x8  mask  (dot- 
dot)  with  the  visible  piece  of  target.  The  4x4  scanning  mask  (dash-dot)  gives 
superior  performance. 
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5.0  FIGURE  OF  MERIT  FOR  THE  CORRELATION  PLANE  PRODUCED 

BY  A  PHASE-ONLY  FILTER 

One  of  the  steps  in  the  Consensus  Correlator  Algorithm  involves  comparing  correlation 
spikes  among  cases  where  different  amounts  of  the  input  scene  are  masked  off  as  well  as  having 
used  different  phase-only  filters.  This  clarifies  the  need  for  a  robust  figure  of  merit.  One  figure  of 
merit  that  was  tested  was  the  ratio  of  the  correlation  spike  height  to  the  secondary  peak  (after 
subtracting  the  background).  This  is  an  important  measure  but  our  implementation  did  not  give 
consistent  results.  A  second  measure  which  showed  some  consistency  is:  (max-mean)/mean.  The 
next  two  subsections  present  our  preliminary  work  on  this  latter  figure  of  merit. 

5.1  Target  Size  Effects  in  Matched  Filter  and  Phase-Only  Filter  Correlators 

To  assess  the  effect  of  the  number  of  pixels  in  a  target  shape  on  matched  filter  and  phase- 
only  filter  correlators,  a  test  using  target  shapes  with  1  to  1024  nonzero  pixels  was  designed  (using 
a  32  x  32  matrix).  The  target  shapes  use  pixels  with  values  zero  or  one,  and  these  are  randomly 
located.  In  this  test  a  new  fourier  plane  filter  that  was  made  to  match  each  random  target  was  then 
correlated  with  its  target.  The  resulting  correlation  peak,  mean,  (peak-mean)  and  the  (peak- 
mean  )/mean  are  graphed  in  Figures  14  and  15  as  a  function  of  the  number  of  pixels  in  the  random 
targets. 

The  Consensus  Correlator  is  designed  to  handle  the  situation  where  smaller  numbers  of 
pixels  are  available,  which  suggests  that  we  examine  the  left  edge  of  the  eight  graphs  of  Figures  14 
and  15  for  constant  behavior.  We  observe  that  (max-mean )/mean  is  the  most  constant,  but  only  for 
phase-only  correlators.  This  result  tentatively  supports  the  (max-mean )/mean  formula  as  a 
reasonable  figure  of  merit. 
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Figure  14.  The  four  graphs  measure  the  performance  of  a  2-d  correlator  where  the  input  is  a 
random  pattern  of  ones  and  zeros  and  the  fourier  plane  Filter  is  the  corresponding 
MATCHED  FILTER.  The  horizontal  axis  is  the  number  of  pixels  that  are  ones. 
The  graphs  beginning  at  the  top  of  the  page  are:  the  correlation  plane  maximum, 
mean,  (max-mean),  and  (max-mean )/mean. 
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Figure  1 5.  The  four  graphs  measure  the  performance  of  a  2-d  correlator  where  the  input  is  a 
random  pattern  of  ones  and  zeros  and  the  fourier  plane  filter  is  the  corresponding 
PHASE-ONLY  FILTER.  The  horizontal  axis  is  the  number  of  pixels  that  are 
one’s.  The  graphs  beginning  at  the  top  of  the  page  are:  the  correlation  plane 
maximum,  mean,  (max-mean),  and  (max-mean )/mean. 


5.2  Correlation  Plane  Statistics  for  Uncorrelated  Noise  Inputs  as  a  Function  of  the  Size  of  the 
Expected  Target 

An  evaluation  of  the  behavior  of  the  measure  (max-mean)/mean  for  real  target  filters  but 
with  pure  unrelated  noise  as  the  input  is  necessary.  This  is  important  because  it  provides  a  more 
absolute  measure  of  the  significance  of  a  given  correlation  spike. 

This  was  investigated  by  using  a  random  noise  input  (absolute  value  of  Gaussian  noise) 
and  a  phase-only  correlator  whose  filter  is  derived  from  a  real  target  shape.  This  real  target  shape 
was  based  on  our  standard  car  target  but  with  a  variable  number  of  pixels,  starting  with  just  one 
column  of  pixels  at  the  left  and  increasing  the  number  of  pixels  by  including  progressively  more 
columns  stepwise  to  the  right. 

Figures  16  and  17  show  the  mean  and  standard  deviation  of  the  measure  [(max- 
mean  )/mean]  as  a  function  of  the  number  of  pixels  in  the  basis  target  of  the  phase-only  filter,  using 
5 1 2  different  input  noise  fields.  The  dashed  lines  are  the  data  and  the  solid  lines  are  a  simple 
exponential  model. 

The  data  in  these  two  subsections  suggests  that  the  measure  (max-mean)/man  is  somewhat 
insensitive  to  the  number  of  pixels  in  the  basis  target  for  target  detection  and  follows  a  predictable 
course  for  noise  inputs.  The  figure  of  merit  tentatively  adopted  for  (max-mean)/mean  was  the 
number  of  standard  deviations  above  the  noise  norm  of  [(max-mean )/mean]. 

This  figure  of  merit  was  used  in  the  algorithm  with  reasonable  success.  However  the 
effects  of  the  thresholding  decisions,  the  larger  number  of  cases  occuring  with  the  4  x  4  pixel  mask 
over  that  of  the  16  x  16  pixel  masks,  and  the  regenerative  optimization  used  on  the  fourier  plane 
filter  have  changed  the  statistics  and  made  the  simple  figure  of  merit  inadequate  for  a  real  system. 
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Figure  16.  Correlating  a  phase-only  filter  which  was  designed  to  detect  a  piece  of  a  known 
target  (scout  car)  with  an  unrelated  noise  field  shows  that  the  512  case  average 
value  (vertical  axis)  of  [(max-mean)/mean]  of  the  correlation  plane  decreases  as  the 
number  of  pixels  (horizontal  axis)  increases  in  the  intended  target  shape. 


Figure  17.  This  gives  the  standard  deviation  (vertical  axis)  of  [(max-mean )/mean]  of  the 

correlation  plane  versus  the  number  of  pixels  in  the  intended  target  which  is  used  to 
create  the  phase-only  filter.  The  input  is  uncorrelated  noise. 
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6.0  A  COMPARISON  OF  MATCHED  FILTERS,  PHASE-ONLY  FILTERS 
AND  HYPER  PHASE-ONLY  FILTERS 


The  unusual  situation  of  our  Consensus  Correlation  where  most  of  the  target  shape  is 
missing  suggests  we  review  the  relative  merits  of  the  phase-only  filter  and  the  standard  matched 
filter. 

6. 1  A  Comparison  of  a  Matched  Filter  and  a  Phase-Onlv  Filter  tor  Detecting  a  Small  Piece 
of  the  Expected  Target 

A  simple  test  was  performed  using  our  standard  target  to  generate  a  matched  filter  and  a 
phase-only  filter.  The  input  image  to  the  correlator  was  a  masked  off  version  of  the  target  shape, 
passing  only  a  4  x  4  oixel  area  of  the  target.  Figure  1 8  shows  the  correlation  plane  using  the 
matched  filter.  There  are  three  peaks  of  equal  height.  This  means  that  the  masked  off  version  of 
the  target  shape  found  3  locations  with  the  original  target  shape  that  gave  100%  correlation. 

Figure  19  shows  the  correlation  plane  of  the  same  masked  target  but  using  a  phase-only 
filter.  It  has  two  “incorrect”  spikes  which  are  twice  as  large  as  the  spike  which  correctly  locates  the 
target.  We  observe  then  that  a  phase  only  filter  can  have  erratic  behavior  for  highly  masked 
targets.  This  is  one  of  the  reasons  that  we  arranged  for  the  phase-only  filter  in  the  Consensus 
Correlator  Algorithm  to  be  recalculated  to  better  match  the  visible  pieces  of  the  target. 

6.2  Spectral  Whitening  of  the  Input  Target  Gives  the  Sharpest  Correlation  Spike 

This  subsection  presents  a  series  of  correlation  planes  obtained  using  the  car  target  with  a 
sequence  of  fourier  plane  filters  ranging  smoothly  from  a  matched  filter  to  a  phase-only  filter  and 
further  into  the  realm  of  hyper  phase-only  filters. 
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Figure  1 8.  The  correlation  plane  for  a  4  x  4  pixel  piece  of  the  scout  car  target  as  the  input  and  a 
matched  filter  based  on  the  whole  target  in  the  fourier  plane.  Three  spikes  of  equal 
amplitude  appear. 


Figure  1 9.  The  correlation  plane  for  the  same  4x4  pixel  piece  of  the  scout  car  used  in  the 

previous  figure.  Here  a  phase-only  filter  based  on  the  whole  target  shape  gives  two 
spikes  twice  as  large  as  the  actual  correct  correlation  spike. 


Denoting  a  complex  array  mf  for  matched  filter,  we  create  a  phase-only  filter :  pof  = 
mf/[e+lmfl]  where  e  prevents  dividing  by  zero  and  we  understand  that  this  is  an -element  by  element 
array  division.  This  can  be  generalized  to  a  general  filter:  gf  =  mf/[e+|mfl]exP.  For  exp  =  0  we 
create  a  matched  filter.  For  exp  =  1  we  create  a  phase-only  filter.  For  exp  =  2  we  create  a  hyper 
phase-only  filter;  which  is  actually  not  a  phase-only  filter  since  it  has  amplitude  information. 

Figures  20  and  21  show  a  sequence  of  correlation  planes  beginning  with  a  matched  filter 
(exp  =  0)  and  smoothly  changing  with  exp  =  0,  .25,  .50,  .75.  1.0,  and  2.0.  The  correlation 
spikes  get  sharper  as  the  exponent  approaches  2.0.  The  exponent  exp  =  2  has  the  effect  of 
whitening  the  spectrum  of  the  target  input,  all  spatial  frequencies  that  are  nonzero  now  have  the 
same  amplitude.  This  extreme  adjustment  should  also  have  remarkably  bad  noise  properties. 
However  in  the  complete  absence  of  noise  the  ratio  of  the  correlation  spike  height  to  the  secondary 
peak  can  achieve  very  large  values.  Figure  22  shows  the  correlation  spike  height  to  secondary 
peak  for  the  car  target  as  a  function  of  the  exponent  discussed  above.  The  matched  filter 
corresponds  to  exp  =  0;  the  phase-only  filter  has  exp  =  1  and  the  fully  whitened  target  spectrum 
has  exp  =  2. 
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Figure  20.  The  top  correlation  plane  is  the  result  of  the  scout  car  as  input  and  a  matched  filter. 

The  matched  filter  has  exp  =  0  in  gf  =  mf/[e  +  |mf  |]exp.  The  middle  correlation 
plane  uses  exp  =  .25,  and  the  lower  correlation  plane  uses  exp  =  .50.  A  phase- 
only  filter  would  have  exp  =  1. 
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Figure  22.  The  ratio  of  the  correlation  spike  height  to  the  secondary  peak  is  shown  versus  exp 
(see  Figures  20  and  21).  At  the  left  with  exp  =  0  is  the  matched  filter,  exp  =  1  is  a 
phase-only  filter,  and  exp  =  2  uses  the  fourier  plane  filter  to  whiten  the  input  target 
spectrum  resulting  in  an  impulse-like  spike  in  the  correlation  plane. 


7.0  CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  successfully  demonstrated  the  Consensus  Correlator  for  detecting  targets  that  are 
mostly  obscured.  The  Consensus  Filter  mimics  the  human  capacity  to  notice  pieces  of  a  target  and 
to  associate  them  together  to  determine  if  they  could  be  part  of  a  single  target  shape,  while  ignoring 
the  interspersed  clutter.  This  demonstration  used  random  noise  as  the  background.  We 
recommend  that  the  second  demonstration  be  on  targets  that  are  mostly  obscured  and  have  a 
camouflage  background  consisting  of  randomly  situated  pieces  of  the  actual  target.  The  Consensus 
Correlator  will  have  the  advantage  here  of  removing  those  camouflage  pieces  which  are  not  situated 
properly  to  compose  a  real  target. 

The  demonstration  involved  a  complicated  computer  program  which  we  expect  would  have 
been  even  more  successful  had  it  been  optimized.  A  portion  of  the  work  focused  on  testing  figures 
of  merit  which  is  a  central  issue  for  complex  optical  correlators.  Our  progress  in  figures  of  merit 
was  useful  to  our  own  Consensus  Correlator  demonstration,  and  additional  work  in  this  area 
would  probably  receive  broad  interest  from  other  researchers  in  this  field. 

Our  demonstration  used  phase-only  filters  which  facilitated  our  finding  a  useful  figure  of 
merit  but  also  produced  anomalies  in  the  detection  of  mostly  obscured  targets.  We  partially 
compensated  for  this  behavior  by  recalculating  the  phase-only  filter  using  estimates  as  to  what  parts 
of  the  target  were  actually  present. 

Another  avenue  for  increasing  our  understanding  of  phase-only  filters  would  be  to  examine 
the  phase-only  filter  in  the  convolution  domain:  that  is,  inverse  fourier  transform  a  variety  of  phase- 
only  filters  to  get  experience  with  their  properties  in  the  image  domain. 
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8.0  APPENDIX 


The  appendix  contains  a  list  with  a  short  description  of  the  important  functions  and  script 
files  in  Section  8. 1 .  Section  8.2  is  a  list  of  the  important  variables.  The  program  codes  developed 
by  the  principal  investigator  are  listed  in  Section  8.3. 

8.1  Functions  and  Script  Files 


run4.m 

scan2.m: 

test3.m 

finmask2.m 

optfilt2.m 


expand(q) 

peaks(x) 

optcor(x,y) 

msk(a,r,c,bh,bw) 

plt4(x) 


Runs  whole  program 

Moves  mask  over  input  scene  and  collects  correlation  spikes 
Groups  correlation  spikes  from  all  the  cases  of  the  masked  input  scene 
Assembles  masks  based  on  groupings  of  correlation  spikes  and  compares 
the  final  performances 

Calculates  a  phase-only  fourier  plane  filter  optimized  for  the  pieces  of  the 
target  shape  which  have  been  found.  Removes  parts  of  mask  which  cannot 
overlap  with  apparent  target  position. 

Expands  mask  q  assuming  32  x  32  matrix  (4x4  becomes  6x6) 

Converts  matrix  x  to  matrix  with  only  local  maxima  (4  nearest  neighbors) 
Gives  correlation  plane  for  input  scene  x  and  conjugate  fourier  plane  filter  y 
Makes  ones  in  matrix  a  from  (rows)  r  to  r+bh- 1  and  columns  from  c  to 
c+bw-1. 

Plots  four  graphs  from  matrix  with  five  columns  with  the  x  coordinate  being 
the  first  column  and  the  other  four  columns  being  the  y’s  of  the  four  graphs. 
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8.2  Variables 

ms 

bh 

bw 

in 

cfni 

thrsl 

spkbx 

ff 

rto2 

rto3 

rfo4 

erne 

esg 

sndpk 

cpl 

space 


matrix  size,  i.e.,  ms  =  32  for  a  32  x  32  matrix 

box  height,  i.e.,  bh  =  4  for  box  height  of  4  pixels  (chops  up  input  scene) 

box  width 

input  scene  matrix 

conjugate  of  the  2-D  fourier  transform  filter 

threshold  1  =  minimum  spike  to  secondary  ratio  for  all  but  largest  spike 
spikes/box  =  maximum  number  of  correlation  plane  spikes  accepted  per 
chunk  of  input  scene 

attenuation  of  acceptable  subsequent  spikes  below  largest  spike 

ratio  of  primary  correlation  spike  to  secondary  spike  (1  st-median )/(2nd- 

median) 

(max-mean)/mean  in  correlation  plane 
(rto3-eme)/esg 

theoretical  (estimated)  mean  value  of  rto3 
theoretical  (estimated)  standard  deviation  of  rto3 
second  peak  value 
correlation  plane  (matrix) 

allowable  separation  of  correlation  spikes  to  be  considered  fitting  together 
into  one  target 
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’•RUN4.  *!  gene'-'at  es  aat  3h  <  bh  >  ns  ( ns  >  =  Crow,  co  1 ,  rt  o,  bxs,  bh,  rt  o3H 

■/■t  1  =c  1  ock  ; 
c.!  at  4=  Z1 
cat  5= Cl 

for  ns=0:10  %gen  dat4  matrix  with  rto£  arid  dat5  with  rto 
i  r  i  =  c  a  r  8  +  ( n  g  *  0 .  £  '  *  c  a  a  s  n  o  i  5  e ;  input  noise  background 

d  at  4  < ns+ 1,1) =ns ; 
daJ:5  (ns+1 ,  1  )  =ns  ; 
b‘-'=c;*riis ; 
for-  d  i  v=  1  :  4 
bh=bh/£ 
ljw=bh  ; 
s.  can8 
Iest3 
*  i  ninaskil' 

;  '  ( abs i row- 1 ~ '  =space  &  ab? (col -17)  <=5cacp - 

dat  4 ( ns+ 1 , d l v+ 1 > -rt o 4 ; 

dat  5  ( ns+  1 ,  d  i  v-*- 1  )  =- 1  o  : 

dat5(ns  +  l,  div  +  1  '  =rnxdvfa  ; 

end 

end 
eri  c:i 

save  cons7.  mat  dat 4  dat5  space  spkbx  f*  cfT:  t-v-sl 
*/-t  c=c  1  ock  ; 
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■/.SCONE.  M 

■/.  Generates  dat  1  using  ff,  ms,  bh,  bw,  in,  cfm,thrsl, 
'At  -  cl  ock ; 

cl  at  1  =  zeros  (  (  f  1  oor  ( ms  /  bh  )  )  *snkbv.  *  <  f  1  oor  ( ms  /  bw  >  )  ,  7 ) 

i  l  =0 ; 

for  rn=  1  :  f  1  oor  ( rns /  bh  )  , 
for  n=l  :  f  1  oor  (rns/ bw)  , 
mk  =  zeros  (rns,  rns  >  ; 

rnk— rnsk  ( rnk ,  1  +  ( rn  —  1  )  *bh ,  1  +  ( ri—  1  )  *bw,  bh  ,  bw )  r, 
col=optcor  (rnk.  *\n,cfni>  ; 
mesh ( cd 1 ) 

me=mean  ( mean  ( cp  1  )  )  ; 
sgm=std (cpl  (  :  )  )  ; 
rnx  x=rnax  < max  (cpl)  )  ; 
rto=  ( rnx x— me )  /sgm  ; 
f f rto=f f*rto ; 
for  22=1 :spkbx, 

[ y y ,  kl  =rna x  (cpl  )  ; 

Crnx,  col3=max  <yy)  ; 

row=k ( co 1 )  ; 

rto=  ( mx-me  )  /sgrn  ; 

l  f  (rto  <  =  f  f  rto)  ,  break,  end  ; 

col  ( row,  co !)=(>; 

i i=i i  +  1 

dat  1  ( l  i  ,  1:5)=  Cm,  r.,  rn>: -me,  row,  co  1  j  ; 

l  f  ( rto  <  =  thrs  1. )  ,  break  ,  end  ; 

end 

end 

end 

rtat1=datl(l:ii.  :  ) 


spk.bx 
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TEST3.  Uses  dot  1  and  space  (0  to  4)  for  initial  "sec" 

•/.This  is  a  script  file  version  of  tests,  m  and  changes  dat  1 
•/.  It  finds  neighbors  for  Minimum  seperat  ion  "sec" 

Cll,  c]=sue(datl)  ; 

for  sp=space:5 
y -zeros (11, 1 ) ; 
dat 1 ( : , 7) = zeros (11, 1 ) ; 
cor  pp=l:ll 

>•:(:,  1  )  =  (abs  (dat  i  (  pp,  4>  -dat  <=sp  &  abs  (dat  1  (  pp,  5 )  -dat  1  . ,  2  : 

.«  ( : ,  1 )  =  ( x ( : ,  1 ) == 1  ®  (dat 1  (pp,  1 ) ~=dst !(:,!>  1  dat 1 ( pp, £ ) ~=dat 1  i  . , 2 

dat  1  C  : ,  7)  =dat  1  (  ; ,  7)  +  >:  (  : ,  1  )  ; 

■/.  bxs— ■  0;  "/this  section  gives  number  of  different  boxes 

•/.  w  =  C-l  —11;  '/-coritairi  neighboring  spikes  for  each  spike 

■/:  for  qq  =  1  :  1  1 

■/.  if  ( x ( q  q ,  1 > -=  1  f  (any  (dat  1  (qq,  1:2)  ~=w>  )  ) 

’/  bxs=bxs+ 1  ; 

%  w=dat 1 (qq, 1:2) ; 

■/.  end 
%  end 

•/.  dat  1  (pp,  6)=bxs; 
end 


sep=sp; 

i f  (max (dat 1  (  : ,  75  ’  )  0)  ,  break, end ; 
end 

Cmx  rowl=max (dat 1  (  :  , 3)  )  ; 
if  dat 1  ( row,  7 ) =  =  0 
dat 1 (row, 7) -1 ; 
end 

for  pp= 1 : 1 1 
i f  dat 1 ( pp, 7) ==0 
dat  1  (  pp,  7 )  =  1  ; 
end 
end 


t  ha 
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*  PINM«SK£. M 

viThis  is  a  script  file  that  assembles  masns  and  compares  them 
dat£=dat 1 (dat 1 ( : , 7) > 0,  :  )  ;  "/.removes  rows  without  neighbors 
%dat£“dat 1 (dat 1 ( : ,  6) > 0,  : )  ;  Sremoves  rows  without  neighbors 
%dat£ (  : ,  7) =dat2 (  : ,  £)  ; 

Cll,  r3=si:e(dat2)  ; 
if  1 1==0 

CrrtxS,  r  1  =maK  (dat  1  (  :  ,  3)  )  ; 
dat2=dat 1  (r,  :  >  ; 

1  1  =  1; 
end 

hh=f loor (ms/bh ) ; 
ww=f  loor  (rns/bw)  ; 
sx=-2*ones < hh*ww, 1) ; 
x2  =  zeros (hh*ww, 2)  ; 
for  m= 1 : hh 
for  n= 1 : ww 

x£ ( ( m— 1 ) *ww+n,  1  :  2 )  =  C m  n 3  ; 
end 
end 
rsp=0; 

Cr.ix7,  r3=max  (dat2  (  : ,  7)  )  ;  y-finds  max  number  of  neighbors  and  the  first  row 

1 1  =0 ; 
rp=0; 
rpo=0; 
rpo3=0; 
r po4=- 1 000 ; 
mxd vf a= 1000*eps ; 

Awhile  (max (dat2 ( : ,  7)  )  > =0.  5*mx7)  ;  %  all  groups  having  half  as  many  neighbors 

while  (max <dat2 ( : ,  7) ) > =1  !  max (dat2  ( : , 7) ) > *0. 5*mx7) ; 

tt=tt+l ; 

x=*zeros  (11,  1 )  ; 

Cmx77, r3  =max (dat2 ( : , 7) )  ; 

x  ( ; , 1 )  =  ( abs ( dat £ (r,  4) — dat 2 ( : , 4) )  <=sep  4  abs (dat  £ ( r,  5> -dat £  ( : ,  5)  )  (  =  seD)  ; 
x  < : , 1 >  =  (x  < : , 1 ) ==1  «  (dat£(r,  l)~=dat£(;, 1)  :  dat£ (r, 2) ~=dat£ (:,£))) ; 

srnx=sum  (  x  (  : ,  1 )  )  ; 
x (r, 1 ) =1 ; 

x£ ( : , 3) =zer os <hh*ww,  1 )  ; 
rnx3d*  (max(dat£(x,3)  )  )  /10; 
for  pp=l : 1 1 

if  (x(pp, 1)==1  S  dat2 ( pp, 3) > =mx3d > 
m=dat£ ( pp, 1 ) ; 
n*dat£  ( pp,  £)  ; 
x£( (m-1 ) *ww+n, 3) =1 ; 
end 
end 


i  f  any  (  ( x2  (  : ,  3)  ’  *sx=>=sum  (sx  )  >  4  ( sum  (sx )  ==sum  <  x£  ( : ,  3)  >  )  ) 
if  smx>0 
x (r,  1 ) =smx ; 
end 

d«t2(r,  7) =dat2 (r,  7) -x  (r,  1 )  ; 

else 

rsp»rsp+l ; 

sx ( s ,  rsp) =  x2 ( : , 3)  ; 

«ik«z»roB  (ms,  ms)  ; 
for  sss*l»hh*ww 
if  xS(sis,3)xxl; 
m«x2 (sss,  1 )  ; 
n*x2(sss,  2)  ; 

nkassk(mk,  1+ (m- 1 > *bh,  1-Mn-l  >  *b*»,  bh,  bw>  ;  *adds  box  to  mask 
end 
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end 


X  contour  (rnk,  1) 
cp  1  =optcor  (mk.  * i n,  cfm  )  ; 

X  mesh  (cpl ) 
if  smx>0 

X  ( t- ,  1  )  =s mx  : 

end 

datc(r,  7)=dat£(r,7)-x(r,  1)  ;  •/.  reduces  neighbor  number 

7.  dat£ 

[yy,  k]  =rnax  (cpl  )  ; 

[mx ,  col  3  =max  < yy '  ; 
row  =  k (col )  ; 
pts=148; 

if  (bh'-=3a  ;  bw~=3a: 

optfilta  V.  optimizes  the  mask  and  filter 

end 

"A  [row ,  co  1 1 

f  a=  ( sum  ( sum  ( ink  )  )  )  /  (ms  3)  ;  'A  mask,  fractional  area  (noise  est.) 

me=rnean  (mean  ( cp  1  )  5  ; 
rto3=  (mx-rne)  /me; 

eme  =  7  1  7.  78*pt  s.  (-0. £3)  ;  •/■expected  value  o c  -'t  o3 
esg=255*pts.  (-0.72)  ;  ‘/.standard  deviation  of  rto3 

rt o4=  ( rt o3-eme ) /esg  ;  V.num  of  std  above  expected  value 

sgm=sqrt  (sum  (sum  (col .  "2)  /  (ms  '£>  )—  me  2  )  ; 

•-  rt o=  (mx-me )  /sgm ; 

•-  :  e  (rto2>r_  5  .-■-/( pks  '  2)  '  =mxdvbx /25>  ; 

X  if  (rto)rpo  S  mx  / ( bxs ' 2) > =mxdvbx /25) ; 

•■  i  f  (rto3)  roo3  $  mx:  1.0e-30)  ; 
if(rto4>rpo4  £  mx > 1 . 0e-30 ) ; 

•/.  if  (mx-me )  /  f  a>  mxdvfa ;  •/.  signal  to  area  (-•oise'  estimate 

cont  our ( mk,  1 ) 

pause (2! 
pk.s=peaks  (  cp  1  )  ; 

*/.  mesh (pks) 
pks (row, col ) =0 ; 
row2=-5 ; 
col2=-5 ; 

for  sk= 1 :  (  1 +2*space )  2 ; 

Cyy2, k2 ] =ma x  (pks)  ; 

Cmx2,  co 121 =max  <  yy2 )  ; 
row2=k£ ( col  2)  ; 

if  (abs  < row-roxc >  >  space )  !  (abs( co 1 -co 12 ) >  space ) , break,  end ; 

pks (row2, col2) =0 ; 

end 

sndpk=max (max ( pks) ) ; 
med=median (median (cpl ) ) ; 
rto2= (mx-med ) / ( sndpk-med ) ; 

%  eval <  C’ dat  3h’ ,  i nt 2st  r ( bh ) , ’ ns’  ,  int£str (ns ) , * =  Crow, col,  rto,  bxs,  bh,  rt  o£3  ’  3 ) 
m*h  (cpl ) 
pause (2) 

Crow, col  3 
rt  o4 

X  dat2(r, 4:5) 

mxdvfa«  (mx-mr*)  /fa ; 
rpo3*rto3; 
rpo4»rto4 ; 
rp-rfc©2j 
rpo“rto ; 
rowlTewj 
col l^col 5 
end  xiE 


rtog»r5y 
rowroBI  < 

c el -coil  • 

rto=rpo ; 

rto3»rpo3 ; 
r t  o4  =  rpo4 ; 
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Opt  f  i  lt£.  rn  is  a  script  file  which  calculates  a  phase-only  £-d 
optical  correlation  filter  which  is  optimized  for 
the  pieces  of  the  target  shape  which  have  been  found. 
Uses  row, col  of  correlation  spike  which  is  expected 
at  17, 17.  Removes  parts  of  mask  that  doesn’t  overlap 
with  detected  target  location. 

rnk£=shift  (rnk,  17-row,  17— col)  ; 
riik£=carc.  *expand  (rnk£)  ; 
cfrii£E:=conj  (fft£(rnk£)  )  ; 
cfr,i£  =  cfnit.  /  (eps+abs  (cfni£)  )  ; 
rn k£— expand  (expand  (rnk£)  )  ; 

cp  1  -opt  cor  (  :i.  n.  ♦rnk.  *sh  i  ft  (mk£,  row- 17,  co  1  -  1 7 )  ,  cfni£)  •; 
r  y  y ,  k 1 =ma  x ( c  p 1 )  ; 

Crux ,  col  1  =max  ( yy )  ; 
row  — U  ( col )  ; 

:ih£=car£.  ♦shift  (rnk,  17-row,  17-col)  ; 
cfni£=conj  (fft£(rnk£)  )  ; 
c'fni£=cfni£.  /  (eps+abs(cfni£)  )  ; 
mk,::'=ex pand  (expand  (rnk£)  )  ; 

Mk-rtik.  *sh  i  ft  (rnk£,  row-17,  col-  17)  ; 
c p  1  =opt cor ( m.  *m« , cf n i £ )  ; 
r y y , k 3 =ma  x  <  c p 1 )  ; 

Crux,  co  1 1  =ma  x  ( y  y )  ; 
row=k v  co 1 )  ; 


■/ 


•/. 

■/. 

•/ 
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function  z=exDand (q) ;  SExpanos  mask  m  32,32  matrix  <4  by  4  into 
a=zeros < 33,  32)  ; 

a=q  +  sh  i  f  t  ( q,  - 1 , 0)  +  sh  i  f  t  ( q,  1 ,  0)  ; 
a=a  +  shi  ft  (a,  O,  — 1 )  ■+•  sh  i  f  t  (a,  0,  1  >  ; 
z=a>  0.5; 


•/-PEAKS  converts  input  matrix  to  matrix  with  only  local  maxima 

"/■  —all  others  set  to  zero,  including  the  outer  rim 

'/•  (maxima  must  be  greater  than  or  equal  4  nearest  neighbors) 

function  y=oeaks(x)  ; 

Cr,  cl =si ze  <  x )  ; 
a=zeros <r, c )  ; 

a  (2  :  r — l,£:c-l)  =  (x<£:» — 1 ,  £ :  c- 1 )  >  =x  (  1 :  i — 2,  2:  c-  1 )  )  ; 
b=a ; 

a  (2  :i — 1, 2 :  c-1 )  =  <  x  <2:  t — 1,  2:  e-1 )  >  =x  (3:  r,  2:  c- 1 )  )  ; 
b=b. *a ; 

a  <  2 :  r — 1 , 2 : c-1 >  =  (  x (2 :  i — 1, 2: c-1) >  =  x (2: r-1 ,  1: c-2) )  ; 

□=b. *a : 

a  (2  jr-l ,  2  :  c-1  )  =  (  x  (2: ) — 1 ,  2 :  c- 1 )  >  =x  (2 :  r- 1 , 3 :  c)  )  ; 
o=o.  a ; 
v  =  d  .  *  x  ; 


function  :=optcor ( x, y ) 

z=  <  abs  <  f f t sh i f t  ( i f ft 2 ( f f 1 2 (x).*(y))))).  ’2 ; 


fc  by  3  > 
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C:\MATLAB\SIGNADtype  rnsk.  m 
function  z=msk  (a,  r,  c,  bh,  bw) 
a  r :  r+bh  —  1 ,  c  :  c+bw- 1  )  —ones  (  bh  ,  Dw )  ; 
z=a; 


function  y=plt4(x)  '/.plots  matrix  as  four  graphs 

"/•assuming  5  columns  with  the 
‘/■X  coordinate  being  the  first 
"/.column  and  the  other  4  columns 
*4being  the  4  graphs  (y’s) 


dat  =x  >  0.  1 ; 

X  (#^cStj)™K  V 

dat  =dat  ==0 ; 
x ( : , 2:5) =x ( 
x 1  =  x  < : ,  1 )  ; 
y  1  =  x  (  :  ,2 )  ; 
y£=x  < : , 3)  ; 
y3=x ( : , 4)  ; 


"/.pad  with  zeros 
2:5).  *dat  < : , £ : 5)  ; 

£ : 5 ) +0.  1 *dat ( : , £ : 5)  ; 


y4=x ( : , 5) ; 

semi  logy (xl,  yl,  xl, y£, 1  — ’ , xl,  y3,  ’  :’ , xl, y4, 


’  ) 


--US  GOVERNMENT  PRINTING  Off  ICE 
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MISSION 

OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 

O 

Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C  I  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


